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Abstract 



_ The amplification of disk oscillations resulting from non-linear resonant couplings 

1^ I between the oscillations and a warp is examined. The disks are geometrically thin 



and general relativistic with non-rotating central objects. By using a Lagrangian for- 
mulation, a general stability criterion is derived. The criterion is applied to horizontal 
^ ■ and vertical resonances of g-mode oscillations and to horizontal resonances of p-mode 

oscillations. The results of analyses show that g-mode oscillations (including p-mode 
O ' oscillations of n = 0) are amplified by horizontal resonances with the warp. Other 

O ! modes of oscillations with other resonances are all damped by resonances. The ampli- 

^ I fied g-mode oscillations are located around the radius 4rg, the radius of the epicyclic 

frequency being maximum, i.e., /tmax- The frequencies of amplified oscillations are 
harmonics of Kma^, and well explain the 2 : 3 pairs of observed QPOs for reasonable 
masses (and spins) of the central objects. 
^ ■ Key words: accretion, accretion disks — black holes — kiloherz quasi-periodic 



oscillations — relativity — warps — X-rays; stars 



1. Introduction 

Kiloherz quasi-periodic oscillations (kHz QPOs) have been observed in many X-ray bi- 
naries. They are observed in both neutron stars and balck holes (for review, e.g., van der Klis 
2000). Recently, they were also observed in X-ray pulsars (Wijnands et al. 2003). One of the 
reasons why much attention is being paid to kHz QPOs is that they give a clue for understand- 
ing the structure of the innermost region of relativistic disks, and provide powerful information 
concerning the masses and spins of the central objects. 

One of characteristics of kHz QPOs in black hole binaries is that they appear in pairs 
in many objects, and their frequency ratios are close to 2 : 3. Paying attention on these 
points, Abramowicz and Kluzhiak (2001) and Kluzhiak and Abramowicz (2001) suggested that a 
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parametric resonance is the cause of these QPOs. Subsequently, they and their group developed 
their idea in many subsequent papers (as recent papers, see, for example, Abramowicz, Kluzniak 
2003; Abramowicz et al. 2004a; Abramowicz et al. 2004b). For a parametric resonance to work, 
disks must be deformed. A well-known example of the excitation of oscillations on deformed 
disks is a precession of disks deformed by a tidal force (Whitehurst 1988; Hirose, Osaki 1990; 
Lubow 1991). Another example is excitation of a spiral pattern in galactic disks deformed 
by ram-pressure (Tosa 1994; Kato, Tosa 1994). In the context of kHz QPOs, Kato (2003) 
considered disks deformed by a warp, and studied non-linear resonant interactions between the 
warp and g-mode oscillations on disks. He showed that some g-mode oscillations are amplified 
on the disks by resonances, and suggested that they would be the cause of QPOs. 

The main purpose of this paper is to re-examine more generally and more carefully the 
amplification of disk oscillations resulting from resonant couplings between a warp and the 
oscillations, correcting some errors in Kato (2003). That is, we elaborate formulations while 
studying the condition of nonlinear resonant amplifications by adopting the Lagrangian for- 
mulation (Lynden-Bell, Ostriker 1967) of hydrodynamical equations. [In Kato (2003) Eulerian 
formulations are adopted.] There are two superior points in Lagrangian formulations. The first 
is that the non- linear coupling terms expressed in Cartesian coordinates are compact. The sec- 
ond is that the linear part of hydrodynamical equations can be expressed in terms of Hermite 
operators (Lynden-Bell, Ostriker 1967). This makes stability analyses perspective, and allows 
us to derive a general stability criterion. 

As a preparation for systematically studying which resonances amplify oscillations, the 
types of resonant couplings between a warp and oscillations were examined in detail in a 
previous paper (Kato 2004, hereafter Paper I). The second purpose of this paper is to show, by 
calculating the growth rates, which resonances really amplify oscillations, and to suggest the 
origin of the observed kHz QPOs. 

2. Basic Non-Linear Equations and a General Stability Criterion 

An unperturbed disk is geometrically thin, non- self gravitating, and relativistic without 
accretion flow. The effects of general relativity are taken into account within the framework of 
a pseudo- Newtonian formulation using the gravitational potential introduced by Paczyhski and 
Wiita (1980), no rotation of the central object being assumed. This means that all calculations 
of non-linear coupling terms are made within the framework of the Newtonian hydrodynamics, 
changing only the Newtonian gravitational potential to that of Paczyhski and Wiita. 

2.1. Non-linear Hydrodynamical Equations of Lagrangian Formulation 

In a Lagrangian formulation, the hydrodynamical equation describing perturbations can 
be written as (Lynden-Bell, Ostriker 1967) 
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^'^-6(-V^--Vp), (1) 



Df^ V p 

where ^ is a displacement associated with perturbations, and D^/Dt is the time derivative 
along an unperturbed flow, Uq, and is related to the Eulerian time derivative, d/dt, by 

In equation (1), 5 represents the Lagrangian variation of the quantities in the subsequent 
parentheses, and ip is the gravitational potential. Other notations in equation (1) have their 
usual meanings. 

The Lagrangian change of a quantity, Q{r,t), i.e., 6Q, is defined by 

6Q = Q[r + ^{r,t)]-Qo{r,t), (3) 

where the subscript to Q represents the unperturbed value. The Eulerian change of Q, i.e., 
Qi, on the other hand, is defined by 

Q, = Q{r,t)-Qo{r,t). (4) 

If Q[r + ^{r,t),t\ is Taylor-expanded around Q{r) up to the second-order terms with respect 
to perturbations, we have, from equations (3) and (4), 

This is a relation between 5Q and Qi, up to the second-order small quantities with respect to 
perturbations. Here and hereafter, subscripts and 1 are used to represent the unperturbed and 
(Eulerian) perturbed quantities, respectively. Furthermore, here and hereafter, the summation 
abbreviation is adopted, using Cartescian coordinates. (In the final calculations of the growth 
rate of the resonance interactions we adopt cylindrical coordinates.) 

Our purpose here is to explicitly write down equation (1) up to the second-order small 
quantities with respect to perturbations. Since the self-gravity of the disk is neglected here, 
i.e., iIji = 0, we have easily 

i(V.,)=&|-(V.)„) + if.f.g;^^(V.)o). (6) 

The second term on the right-hand side of equation (6) represents non-linear terms. Expressing 
(5(Vp/ p) in terms of Lagrangian quantities is, on the other hand, somewhat complicated. Using 
equation (5) and the definition of Lagrangian change, we have 

d_ 

The Eulerian pressure variation, pi, which appears on the right-hand side of equation (7), is 
expressed in terms of 6p and $, by using 
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Vpo + Vpi + 0— V(Po + Pi] 



6p - 

We then have 



dridvj 



(8) 
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5p\ 
Po/ 



(9) 



This exphcitly expresses p^diVp/ p) up to the second-order quantities in terms of Lagrangian 
qunatities, 5p, 6p, and ^. 

Next, 6p and 5p on the right-hand side of equation (9) are exphcitly expressed in terms 
of ^. To do so we use the equation of continuity and the adiabatic relation. The former relation 
is expressed as (e.g., Kato, Unno 1967) 



5p + pG 



dr-i 



PoNc, 



where 



The adiabatic relation is written as (e.g., Kato, Unno 1967) 



Po Vo 



where 
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(11) 



(12) 



(13) 



Substitution of equations (10) and (12) into equation (9) finally gives an expression for 
Po<^(Vp/p), expressed in terms of $, alone. This is summarized as 



Pod ( - Vp 



(1) 



(2), 



(14) 



where is the linear part with respect to $, and P(2) is the non-linear part. An explicit form 
of the linear part is given by Lynden-Bell and Ostriker (1967), which is 



PW = V 



:i-ri)po 



PoV 



dpo 
dr,; 



(15) 



The non-linear term, P{2)i generally has a complicated form. In the case of isothermal disks, 
i.e., Fi = 1, however, it has a simple expression. We are satisfied here with using this simple 
expression for P(2) in the following analyses of resonances, since the difference of Fi from unity 
will not bring about any qualitative differences in the final results. In the case of Fi = 1, after 
some manipulations we have 
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Pi^-^ = ^U^^^i,Y (16) 



Hereafter, we consider time-periodic perturbations, which vary as exp(zu;t), uj being a 
frequency. Then, the non-hnear hydrodynamical equation for adiabatic perturbations, equation 
(1), is written as 

-^Vol + 2z^Po(^^o ■ V)^ + L{i) = poC(^,^), (17) 

where X is a hnear operator with respect to ^, and is (Lynden-Bell, Ostriker 1967) 

Hi) = Po(i^o ■ V)(uo ■ V)^ + po(^ ■ V)(V^o) + P(i)(^). (18) 

The non-hnear term p^C consists of two parts, i.e., p^C = poC^ + poCp, where poC^ comes 
from —pq6(ViP) and poCp from — Po<^(Vp/p). From equation (6) we have 

PoC^ = -^Poe.e,^^(v^o), (19) 

and from equation (16) we have 

PoC. = -^U|^VOy (20) 



Detailed expressions for p^C^ and p^Cp in cyhndrical coordinates are given in appendix 1. 

An important characteristics of the hnear operator L is that it is a Hermitian (Lynden- 
Bell, Ostriker 1967). This characteristics of L is used later when we evaluate the growth rate 
of oscillations. 

2.2. A General Method for Deriving Stability Criterion 

We assume that a mode of oscillations is present on the disk with a small amplitude. The 
displacement associated with the oscillations is denoted by ^, which satisfies the homogeneous 
part of the wave equation (17). Now, we further assume that the unperturbed disk is deformed 
by some external forces from an equilibrium state. The displacement associated with the defor- 
mations is denoted by Nonlinear couplings between ^ and ^™ are now examined, assuming 
that they are weak and their effects can be studied by applying successive approximations. 
The processes of the couplings are shown schematically in figure 1. The first step of nonlinear 
couplings between ^ and is their direct interactions. The resulting oscillations are denoted 
by and hereafter called intermediate oscillations. Since we are interested here in the case 
when the time scale associated with the deformation is sufficiently long compared with l/cj, 
the intermediate oscillations are described by using the inhomogeneous wave equation (17) as 

-uj^p,t' + 2z^Po(^^o ■ V)r' + L{t') = PoC(^,r)- (21) 

As the next step of nonlinear interactions, we consider interactions between the interme- 
diate oscillations and the deformations This interaction is interesting, since it generally 
feedbacks to the original oscillations of ^. Unless a resonance phenomenon is involved in the 
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Fig. 1. Feedback processes of nonlinear resonant interactions acting on oscillations. The original oscilla- 
tions are characterized by uj, m, and n. Since the warp corresponds to a wave mode of = 0, m = 1, and 
n = l, the nonlinear interaction between them brings about intermediate modes of oscillations of oj, to± 1, 
and nil. To these intermediate oscillations, the disk resonantly responds at a certain radius. Then, the 
intermediate oscillations feedback to the original oscillations after the resonance. This feedback process 
amplifies or dampens the original oscillations, since resonances are involved in the interaction processes. 
(After Paper I) 

feedback processes, the effects of the feedbacks are only to change the frequencies and func- 
tional forms of the original oscillations. If a resonance process is involved, however, there is 
a possibility of the oscillations being amplified (or damped). The Lagrangian formulation of 
hydrodynamical equations by Lynden-Bell and Ostriker (1967) makes this analysis easy. The 
equation describing these feedbacks is given by 

-u^p^i + 2iup^{u, ■ V)^ + L{$,) = poCit\n- (22) 

That is, the left-hand side of the equation is the same as equation (17), while the coupling term 
is PoC'(^"^';^™)- Let us integrate the equation over a volume with relevant boundary conditions 
after multiplying by where the asterisk denotes the complex conjugate. 

Since the operators on the left-hand side of equation (22) are Hermitians, the integrations 

iJpoC{uo-V)^dV and ji*L{i)dV (23) 

are real. The integration of the right-hand side of the equation, on the other hand, has an 
imaginary part when resonances are involved in the interactions. This implies that uj must be 
complex. If we write uj = u-^ + ioJi, the imaginary part of equation (22) becomes 

-2iuj, J ePo[uJ-tiuo-VMdV = Ql poeCiC\ndV. (24) 

Hereafter, is denoted by u without any confusion. From this equation we can calculate 
In the subsequent sections we examine this issue for the cases where the disk deformations are 
due to warps. 
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Equation (24) has a clear physical meaning. The integral on the left-hand side is the 
wave energy, E, times 2/uj [see equation (92) of Kato (2001)]. Hence, the left-hand side of 
equation (24) represents the growing rate of the wave energy, —2iuJiE, times 2/uj. This implies 
that the right-hand side of equation (24) must be the work, W, done on the original oscillations 
by the force C(^™*,^™) during the time interval 2/uj. In other words, the work done on the 
oscillations by resonances, W, during a unit time is 

W='^qJ eCie'^'^OdV. (25) 

This is understandable, since if the velocity associated with the oscillations and the force, C, 
are in phase as an idealized situation, the work done on the oscillations during which the 
displacement becomes ^ (the quarter of one cycle) is ^3 / ^*C(^™*, ^™)(iV. Hence, the rate of 
work (per unit time) is, on average, / $,*C{^"^^ ,^'")dV times 2uj/'k. It is finally noted that the 
wave enegy, E., is expressed later in a more conventional form. 

3. Warps, Oscillations, and Dispersion Relation 

The unperturbed disk is assumed to have no motion except for rotation. The rotation 
is further assumed to be cylindrical. That is, with cylindrical coordinates (r,(y9,2;) whose 2;-axis 
is perpendicular to the disk plane and whose origin is at the disk center, the angular velocity 
of rotation, fi, is represented as = r2(r). When we need numerical figures, we adopt the 
Keplerian form for fi, i.e., = ^ySj'\ since geometrically thin disks are considered here. 

The unperturbed disk is further assumed to be isothermal in the vertical direction. In 
vertically isothermal disks the density is stratified as (e.g., Kato et al. 1998) 

Po(r,2;) =poo('^)exp(^-^^^^, (26) 

where poo is the density on the equatorial plane and the half-thickness of the disk, H ^ is related 
to the angular velocity of the Keplerian rotation, fix, by 

VLl{r)H\r) = ''-^ = cl (27) 
Po 

where Cs is the isothermal speed of sound. 

3.1. Oscillations Modes and Warps 

Oscillations on disks are generally a combination of motions on the equatorial plane and 
those in the vertical direction. These motions are coupled in complicated ways, and cannot be 
treated separately. In the case of the vertically isothermal disks, however, we can approximately 
separate them (Okazaki el al. 1987; see also Kato et al. 1998). In order to make analyses simple, 
we introduce here an approximation. In this case, the z-dependence of the eigen-functions of 
disk oscillations are described by Hermite polynomials, T-Cn{z/H), n specifying the node number 
in the vertical direction. The horizontal displacements associated with the oscillations, ^ are 
then described by 
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and that in the vertical direction, is 

=7^„_iexp[z(^t-m^)]e>)(r). 



(29) 



(28) 



Here, is an abbreviation of the horizontal displacement vector {^r,C,ip) and m is the number 
of arms in the azimuthal direction. The oscillations are now classified by a set of m and n. 
Furthermore, for a set of (m, n), we have two kinds of oscillations, which are p-mode oscillations 
and g-mode ones. The formers are oscillations of {u — rriQ)^ > nQ and the latters are those of 
(a; — mf2)^ < (see e.g., Kato et al. 1998; Kato 2001). The p-mode oscillations start from 
n = 0, while the g-mode ones start from n = l. In the framework of the present approximations, 
the p-mode oscillations of = are purely horizontal. In order to make it clear, the coefficient 
n has been attached to equation (29). 

In addition to the oscillations mentioend above, the unperturbed disk is assumed to be 
deformed by a global standing {uj = 0) pattern of a warp, which can be approximated in 
the lowest approximation to be (e.g., Kato 2003) 



The present warp corresponds to a radially long wavelength perturbation with n = l and m = 1. 
Here and hereafter, the superscript w is used to represent the quantities associated with the 
warp. 

3.2. Dispersion Relation 

As a preparation to solve the inhomogeneous wave equation (17), we discuss here the 
homogeneous part of the equation. The homogeneous part is generally a set of partial differential 
equations. In some cases, however, the homogeneous equation becomes a set of algebraic 
equations, and thus the inhomogeneous non-linear equation (17) becomes manageable. As such 
a case, we discuss the case where the wavelength of perturbations is sufficiently short compared 
with the radius of the disk. This approximation is allowed, except for some special situations, 
since we consider geometrically thin disks. 

Here, we assume that $, oc exp{ikr) and /cr ^ 1. Since the z-dependences of ^ are given 
by equations (28) and (29), using equations (15) and (25) and the characteristics of the Hermite 
polynomials we have, after some manupulations, that the homogeneous part of equation (17) 
is expressed as 



(30) 



and 




(31) 




mfi)2 + K2-4fi2]^"(°)-z2fi(, 



LU — ml 



+ {kH)nl {kH)ii'^ + zef = 0,(32) 




(33) 
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and 



[—{uj — m\ 



(34) 



In these equations the terms with Cg come from the Eulerian pressure variations (see the next 
paragraph). If these terms are neglected, the above equations are those derived by the particle 
approximations, and thus the horizontal motions and the vertical motions are uncoupled. This 
is also shown clearly by deriving the dispersion relation (i.e., the solvability condition) from the 
set of equations (32)-(34), which is 



It will be instructive to note here that for a polytropic gas we can derive after manipu- 
lations that -Pi(^) is related to poV/ii by 



where hi=pi/pQ is the Eulerian pressure variation, pi, being normalized by po- In the framework 
of the present approximations, dhi/dr = cl[k'^^^^'^ + i{k / H)^f^]'Hn and dhi/dz = [nVt\S,f'^ — 
incl{k/H)i^^)\Hn-i. 

4. The First Stage of Coupling between Warp and Oscillations 

We here examine how the disk oscillations and the warp described in section 3 are 
coupled. The procedures of examinations are the same as those in Paper I and outlined in 
subsection 2.2. That is, we assume that the amplitude of oscillations is so small that their 
non-linear effects on the oscillations, themselves, are negligible. The amplitude of the warp, on 
the other hand, is so large that the effects of the warp on the oscillations are non-negligible. We 
examine them successively, as is shown in figure 1. The first step of their interactions is direct 
non-linear coupling between the warp and the oscillations, which brings about intermediate 
oscillations of (mil, 1) with frequency uj (see figure 1). To these intermediate oscillations 
the disk responds resonantly at a radius where the dispersion relation to the intermediate 
oscillations is satisfied. After having this resonant interaction, the intermediate oscillations 
interact with the warp to feedback to the original oscillations of (m, n) with uj (see figure 1). 
This non-linear feedback to the original oscillations through a resonant coupling amplifies or 
dampens the original oscillations. 

In this section we consider the first stage of the above couplings. The intermediate waves 
resulting from this coupling are expressed by attaching a superscript (1) as As mentioned 
before, the coupling has two paths. One is the path where the azimuthal wavenumber of the 
intermediate oscillations is m + 1, and the other is the path with m — 1. To distinguish in 
these two paths, we attach the subscript + to as in the former case and the subscript 
— in the latter case as . From equation (17) the inhomogeneous wave equation describing 
is found to be 



(35) 




(36) 
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-a;Vod'^ + 2iujpo{uo ■ V)d'^ + = PoCT, (37) 

where 

cV'^ = C(£(°\r) and cL'^ = C(^(°\r*)- (38) 

Equation (37) represents two equations. One is the equation for the upper sign (it is + in the 
present case) and the other is that for the lower sign (e.g., — ). Both equations are written in 
one set, since they have similar forms and can be treated in parallel. 

Our purpose here is to solve this inhomogeneous wave equation (37) to obtain 
and Although this equation can be formally solved by using a Green function, it is too 
complicated and not practical. Hence, we are satisfied here to solve the equation by introducing 
the approximations of vertically isothermal disks (see section 3). In this case, concerning the 
functional forms in the vertical direction, the feedbacks to the original oscillations occur through 
two paths. One is a path through n — n — 1 — > n (lower path). The other is through the path of 
n —* n + 1 ^ n (upper path). In the former case, the intermediate oscillations in the horizontal 
direction, $,^]_\ are proportional to T-Cn-i{z/H), while it is proportional to Hn+iiz/H) in the 
latter case. [It is noted that the vertical displacemtnt, associated with the intermediate 
oscillations are proportional to T-Cn-2{z/H) and T-Cn{z/H).] In order to treat these two cases 
separately, we attach the subscript n to as where n represents n — 1 or n + 1. Then, 
the r-, (/?-, and 2;-components of are written as 

= nnexp{t[ujt - (m ± l)^]}eS,^,(r), (39) 

= nnexp{t[ujt - (m ± m}ii!±A^). (40) 

eS,^ = nn-iexp{t[ujt - (m ± m}i^J:lr,{r). (41) 

It should be noted here that in equation (41) subscript n is attached to ^i]±, although it is 
proportional to 7in-i- We adopt this notation in order to clarify that ^r,±,n' '^il^±,n' and C^^^i.n 
are one of the sets of solutions for the same intermediate oscillations. 

Using equations (39)-(41) and the procedures used to derive the wave equations for 
oscillations [equations (32)-(34)], we can write the equations corresponding to equation 
(21) as 

{-[c. - (m ± l)fi]2+(«:^ - 4fi2)}gl, - i2^l[u - (m ± 1)^]^^^^,, 

-(1) -"(1) " (1) 

+ k C^^r,±,h + '>'^^^z,±,h = ^r,±,n5 (42) 

-[uj-im± l)fi]2e^^L,n + -im± imlXn = C'Si.n, (43) 

{-[u - (m ± + nl^iieS,, - zA;n|eS,. = ^'S,., (44) 

where C'r,±,ri and C*<^^±,^i are the parts of Cl^±^f^ and C^^j_^ that are proportional to exp{i{ujt — 
{m±l)ip)}T-Cn{z/H), respectively, i.e., 
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c. 



(1) 



r,±,n 



nn{z/H)fii^Y>{i{^t - (m ± 1)<^)}C', 



(1) 



Similarly, 



(1) 



(45) 
(46) 

(47) 



Detailed expressions for C*r,±,n! C*<^^j_ j^, and C'i^i^ji are given in appendix 2. 
Equation (42)-(44) are easily solved with respect to ^ "* to obtain 



-c 



(1) 

r,±,n' 



(1) 



-c 



±,n 



(1) 



and ^l^L 



(48) 



where 



e2,=-[(^-(m±l)f])^-nfil] 



(1) 

r,±,n 



2fi 



(1) . 



mfi^ <! kH 



(1) 

r,±,n 



2^ Ml) 



[cu-(m±l)fi]2-«:2^(i) 



[cj-(m±l)fi]^ 



2,±,n 



(49) 
(50) 

.(51) 



Furthermore, D± n is defined as 

D±^f, =|[cj - (m ± 1)1]]2 - k2||[cj - {m±l)Q]^ - nfi^j - A;2c2[cu - (m ± l)!]]^. (52) 

An important characteristics of the intermediate mode, is that it becomes infinite at 
the radius where -D±,n = 0. This implies that the disk responds resonantly to the intermediate 
oscillations at these radii. Since the pressure terms in equation (52) are usually small, D± n = 



is realized at radii close to [u — {m± l)f^]^ — ~ or [u — {m± 1)^^]^ — nQ']_ ~ 0. The 



former resonances are hereafter called horizontal resonances, while the latter ones are vertical 
resonances. 



5. Resonant Radii and Resonant Frequencies 

The resonant radii where the disk responds resonantly to the intermediate oscillations 
must be within the propagation region of the original oscillations. This condition has been 
extensively examined by Paper I. We mention here a part of the results. 

There are two kinds of oscillations on disks, i.e., the p-mode and g-mode oscillations. 
Furthermore, resonances have two types, i.e., horizontal and vertical resonances. Combinations 
of them give us three types of resonances: i) horizontal resonances of the g-mode oscillations, 
ii) vertical resonances of the g-mode oscillations, and iii) horizontal resonances of the p-mode 
oscillations. It should be noted that vertical resonances of the p-mode oscillations are absent. 
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As a demonstration, we present here where horizontal resonances of the g-mode oscilla- 
tions occur and how much their frequencies are when n = 1 and the path of the coupling is n 
— ^ n + 1 and n (i.e., n = 2). As equation (52) shows, the horizontal resonances occur near to 

— (m ± 1)^2 ~ K or — (m ± 1)^2 ~ — k. (53) 

The propagation of the g-mode oscillations with frequency u, on the other hand, is bound in the 
region of —n < u — mQ < n [see dispersion relation (35)]. Since we are interested in oscillations 
whose frequencies are near to the boundaries of the propagation region, we require that 

cu — mf2~— K or cj — mf2~K. (54) 

Requirements (53) and (54) are simultaneously satisfied when Vt ~ 2n. That is, in the case of 
the plus sign of m ± 1, they are satisfied when uj — {m+ l)Vt ~ — k and uj — rriQ ~ k, and in the 
case of the minus sign of m ± 1, u — {m — 1)Q ~ n and u — mQ ~ —n. The requirement Q = 2k 
is realized at r ~ A.Ovg when the Schwarzschild metric is adopted, where Vg is the Schwarzschild 
radius, defined by Vg = 2GM/c'^, M being the mass of the central object. The frequencies of 
the oscillations at the radius are rriQ =p n. In the case of m = 2 they are 0.1875 and 0.313 in 
units of (GM/r3)i/2 ^^^^ p^^^^ 

The above discussions are only for the case of the horizontal resonances of the g-mode 
oscillations with n=l and n = 2. Discussions for other types of resonances (Paper I) show that 
resonances can occur at, and only at, the radii where 

K = o, V2n^-n = K, n = 2K, V3n±-n = K, (55) 

depending on the types of the resonances and modes of oscillations. The radii satisfying the 
above requirements in equation (55) are, in turn, 

r = 3rg, r = 3.62rg, r = i.Ovg, and r = GAGrg. (56) 

Since more detailed discussions concerning the resonant radii and the resonant frequencies 
at these radii are given by Paper I, we proceed here to the issue which resonances lead to 
amplifications of the oscillations. 

To proceed further, however, there is an important issue to be noted here, which is the 
sign of the ratio of u; — (m ± l)f2 and u — mQ at the resonance point. This is important to 
discuss whether the resonances lead to amplification or to damping, as will become clear in 
section 7. The sign is negative in the case demonstrated above. Careful examinations show 
that the ratio is always negative in horizontal resonances of all g-mode oscillations, while it is 
positive in vertical resonances of all g-mode oscillations and horizontal resonances of p-mode 
oscillations (except for n = 0). 
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6. The Second Stage of Coupling 



The next issue is to study whether the resonances discussed in the previous sections 
act so as to amplify the original oscillations or not. For this purpose, the non-linear coupling 
between the intermediate oscillations and the warp ^™ is discussed. The coupling brings 
about oscillations of the same form as the original oscillations, i.e., the r- and (f- components 
of the displacements of the resulting oscillations are proportional to exp[i{ujt — ■mif)]T-Cn{z / H) 
and the z-component is proportional to exp[i{ut — nnp)]T-Cn-i{z / H) . This feedback acts so as to 
amplify or dampen the original oscillations, depending on the phase relations at the resonance 
points. 

Now, we start from equation (24). The equation has a clear physical meaning. The 
left-hand side of the equation can be reduced to 

jiu- mn)po{e^*e^ + e^*e^)dv, (57) 

when an approximate relation between ^^P^ and i.e., 

z(cu-mfi)ef +2fie^°^ =0, (58) 

which is the (yj-component of the equation of motion, is used. If vertical integration is performed, 
equation (57) can be written as 

-2iUi{27iY^'H Jiu- ml])pooK?^*e?^ + {n- 1)^^ ]27rrc/r. (59) 

We now remember that the wave energy is known to be defined by (e.g., Kato 2001) 

[In Kato (2001) E is written by using the Eulerian velocity perturbations, Uj. and u^, instead 
of the Lagrangian displacements.] Hence, equation (57) is 

-Ai^E, (61) 

and is related to the rate of change of the wave energy. 

The right-hand side of equation (24) represents the work done on the original oscillations 
through non-linear couplings with a warp. As mentioned before, there are four paths of the 
feedback to the original oscillations, i.e., two paths ofm— *m±l — >m and two paths of n — >■ 
n±l ^ n. We treat these four paths separately. Equation (24) shows that for the coupling of 
m ^ m+1 ^ the rate of energy change of the oscillations is given by 

-^^-E = P,e^*C%{C,%.indV, (62) 



while 



Ai'^E = 53 / Poe^*C%{C,%,ndV, (63) 
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when the path is m ^ m — 1 ^ m. Here, CYfu for example, represents the couphng between 
and ^™*. To make clear that this is the second stage of the coupling and that the path is 
m — > m + 1 — > m, the superscript (2) and the subscript + have been attached to C. 

When integrations on the right-hand side of equations (62) and (63) are performed, we 
should first remember that among various terms of C^'^\ we need only those terms proportional 
to Hn{z/H) for C^f^ and to Hn-i{z/H) for Cp-*. This is because ^^^^ and .^^"^ are proportional 
to l-in{z/H) and T-Cn-i{z/H), respectively, and Hermite polynomials of different orders are 
orthogonal to each other. Second, we should consider that ^'-^-''s have poles at the radii where 
D = is realized, and the imaginary parts of the integrations (62) and (63) come from integration 
around the poles. Related to this, when we calculate C^'^\^^^\ ^'"*) , the term 1/D in ^^^^ can 
be moved outside of C^^^ as {1/D)C^^\C^^\C*)- 

Based on these considerations, we write the C^'^^^s as 



C 



(2) 



,(2) 



a 



(2) 

z,+,h 



1 



D 



+ ,n 



-Hniz/ H)exp{i{ujt - m(^)}C^_+^^(C+,^i, r]* 
-nn-i{z/H)exp{i{ujt - mip}ci%^^{C+,n, V* 



(64) 
(65) 

dr \ 

- (2) 

where the subscript c denotes the value at the resonance radius, and the expressions for Cj^ _,_ ^ 
and Ci%n given in appendix 3. The right-hand side of equation (63) has the same form as 



Then, the right-hand side of equation (62), for example, can be written as 



00 



equation (66), except that C]_^ f^{C,+ h^fl*) 



and CfX^n{C+,n^Vi*) are replaced by C|^/„,^i(C-,n>'^) 
n(C-,n)'^)- Furthermore, Z)+ in equation (66) is replaced by D^. Expressions for 
are also given in appendix 3. 
To perform integration around the poles, we must be careful about the integration path, 
so that causality is satisfied (e.g., Kato 2003). Let us consider the horizontal and vertical 
resonances separately. For the horizontal resonances we have [u; — (m ± l)f2]^ — ~ and 
around the resonance radius, Tc, the first term of the Taylor expansion of -D±,n gives 



(2) 



and C. 

,(2) 



(2) 



C'_L,_,ri, and C^_^n 



r — Tc 



{m±l)n 



where 



dXnVL 



d\nr dlnr J c 

Considering these expressions and the requirement of causality, we have 

1 . . Tr 



D 



-dr ■■ 



±,n 



2{nn\-K^)^\Gn,, 



j-sign[a; — (m ± \)VL\^ 



(67) 



(68) 



(69) 



Here and hereafter, subscript c denotes the values at the resonance radius. 
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In the case of vertical resonances, [u — {m± 1)^^] — nilj_ ~ and around the resonance 



radius, rc, we have 



where 



r — Tc cj— (m±l)r2 

— r 



Gv,±,n=((m±l)fi[a;-(m±l)fi]^ + 2nffix^)j • (71) 
I dlnr V diwr / ) c 



(70) 



Similarly, we have 
1 



D^/' = ^" 2(nl]i-.^^).|Gv,^| ^^g"[" - ^ '^^l- (^2) 
Taking these considerations into account, we can write equation (62) as 

4(-«:2 + nl]2)J^^.+."|c {uj-mn\ E' ^ ^ 

where W is related to the work done on oscillations, and is 

W = n\ ifX^^CfX^^ + n\ eSUC'S,^ + {n-l)\ eT^C'iS,^ (74) 

and 



In the case of m — > m — 1 ^ m, an expression for the growth rate is similar to equation (73). 
The differences are that the arguments of c'^J^ are now C^Jf^ and r). 

Next, we should mention the case of vertical resonances. Derivations of expressions for 
the growth rate are quite parallel with the case of horizontal resonances, and we have 

o _ ^ \n I sign[a;- (m+l)l]]c Wc . . 

We should especially notice that the sign is different from equation (73). 

7. Growth Rate of Resonant Oscillations 

In this section, the growth (or damping) rates of resonant oscillations are calculated 
using equations (73) and (76). As discussed in section 5, there are three kinds of resonances: 
i) horizontal and ii) vertical resonances of g-mode oscillations and iii) horizontal resonances 
of p-mode oscillations. These resonances are discussed here separately. It is noted that the 
oscillation modes of = are p-modes alone. These p-modes are similar in dynamical behaviors 
with g-mode oscillations of > 1. Hence, the p-modes of n = are treated hereafter together 
with g-mode oscillations of > 1. In other words, when we mention p-mode oscillations, they 
are p-mode oscillations of n > 1. Two basic approximations are introduced here. The first one is 
that Cg /c^ is much smaller than the square of the rotational frequency, which means that we 
adopt [kHY ^ 1- This is a natural approximation in treating disk oscillations in geometrically 
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thin disks. Introduction of this approximation imphes that we are interested in resonances that 
occur near to both boundaries of the propagation and resonance regions. That is, we consider 
resonances where [{u — mfi)^ — k^] [{uj — mQ)'^ — nQ\] ~ and {[cu — (m± — K^jKcu — (m± 
— nfl\} ~ are simultaneously realized. The second approximation we adopt here is a 
local approximation of (/cr)^ ^ 1. This approximation is relevant, except for just the boundaries 
of the resonance and progation regions, since we are treating geometrically thin disks. 

7.1. G-Mode Oscillations 

First, we summarize the relations among and in g-mode oscillations. In g- 

mode oscillations, {uj — mVLY — K^ is close to zero, but can be approximated as (k^ — 

nfi^^)] by using the dispersion relation (35). Hence, from equation (34) we have a relation 
between and if^ as 



(0) 



-i{kH)- 



nVL' 



(0) 



A relation between and ^l^^ is, from equation (33), 



(0) 



2Q 



(0) 



(77) 



(78) 



LU — mVt 

These relations show that in g-mode oscillations, ^^^^ is smaller than and ^^'^ by a factor 
of {kH). 



n — 1 — > n (i.e., n = n — 1), 

expressions for ''s are different, we treat these two cases separately. In the case of the upper 



Since in the paths oin^n + l—*n (i.e., n = n + l) and n 



path (n = 77, + 1), negelcting small-order quantities, we have from appendix 2 



Ml) 



i{kH)Qii 



2 m 



2Vt 



V 



c. 



(1) 



-,n+l 



2 m 



dr 



rj + n 



mil r 
drj 

dr oj 



2n 



mf2 r 



(79) 
(80) 



Here, we do not present Cj^^_^_^n+lJ since it is smaller than Crl+,n+i unnecessary hereafter. 
Although C^ll^n+i is smaller than cl^]_^n+i by a factor of kH, it contributes to the final results. 

In the case of the lower path, on the other hand, using expressions in appendix 2 and 
approximations of kH <^ 1 and kr ^ 1, we have 



a 



,n-l 



c, 



(1) 

ip,+,n— 1 



(0) 

z 



c/lnfi2 



dr 



-rj + {n — 1 



drj 
dr 



nikH)e^ + t{n-l)i 



(0) 



The term C^^]_ is not presented here since it is smaller than C^. 
of H/r, and is unnecessary hereafter. 



(1) 



_i and Cj^^^ 



(81) 
(82) 

_i by a factor 
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7.1.1. Horizontal resonances of g-mode oscillations 

To proceed further, we separately consider two coupling paths ofn—^n + l—^n and n 
— ^ n — 1 — ^ n. 

(a) Coupling through the upper path {h = n + 1) 

In this case, we can express C+ n+i terms of by using equations (49)-(51) under 

the help of equations (79) and (80). The results are 

ei„«^.(..).ief{^-K^(§.^:^)} 



2Q 



fi2 



(1) 

z,+,n+l 



in{kH)- 



-k2 + nfi2 ^-.+."+1- 



(1) 



(83) 
(84) 
(85) 



This shows that C^^]. is smaller that c\_\ n+i by a factor of kH, but it is necessary in deriving 
W. 

" (2) 

The next step is to calculate C_^n_^_i by using equations (83)-(85). By consulting the 

expressions for C{^,^) given in appendix 2, we can derive C{^^ \ 'n*)- Considering the order 
of each term, we have 

dr]*" 

-h \fr, ( • : , -, I n -+- r, 

dr 



A(2) - -,f,(kmn'^ t^^^ ^^oV^^) 



- — 7] +n- 

dr dr 



C. 



(2) 



(86) 
(87) 



' (2) 

Here, Cr|, is omitted, since it does not contribute to W . Using the above expressions, we 
finally have, after lengthy manipulations, a simple expression: 



(2) 



+1 



I 

-/«2 + nfi2 IW,n+ll 



(88) 



(b) Coupling through the lower path {n = n — l) 

In this case, C^ + ^.i and C*<[,"|+^„_i are given by equations (81) and (82). Hence, in the 
lowest order of approximations we have, from equations (49) and (50), 



^r',+,n— 1 



-in{kH)Vt\ 
r d\M 



X 



dr 

2n 



r] + 



dr] 



u — {m + l)f2 

(1) 



;(i) 



2n 



V 



dr LU — {m+ l)f2 r 



(89) 
(90) 



In deriving W, we do not need Q + n-i- 



(2) 



Using the above expressions, we can obtain C_,_„_]^. The results are summarized as 
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The (/?-component, C*^ :,_ is not presented here, since it is smaller than Cl^l^^_i and does not 
contribute to W. Taking into account these results, we finally have, after some manipulations, 

Tj/_„lf(0)* A(2) ,^„_lA|t(0)* ^(2) |2 /qoN 

VV —n\ C;,r,+,n-l^r,+,n-l + yji' ^z,+,n-l^ z,+,n-l — _|_ ^(]2 \Sr,+,n-l\ ■ {^■^) 

7.1.2. Vertical resonances of g-mode oscillations 
(a) Coupling through the upper path (n = n + 1) 

Expressions for (7^^^ and Ci^]_ n+i are the same as those for the horizontal resonances 
given by equations (79) and (80). In the present case of vertical resonances, however, expressions 
for the lowest order of approximations of C+L+i are different from those in the horizontal 
resonances. In the expression for Cr + n [equation (49)], the terms coming from C*r,+,n+i and 
C'l^^+^n+i are smaller than the term of C*2^| „_,_;^, since [uj — {m+ 1)^^]^ — n^^ is now small by the 
resonance being vertical. The results of calculations are 



z,+,n+l 

2 

4 f(0) 



ri + n\— + 



dr \dr LO — r 



(94) 



where the expression for (j^^^^^+i is not presented, since it is the same as equation (84). These 
equations show that, different from the case of horizontal resonances, the intermediate oscilla- 
tions have the largest amplitude in the vertical direction. 

By using C+n+i given above, we can calculate c\_\_^_i. The results show that the main 
terms of C^,U, are 

C;5.«=niCT',.«(-^.r + «^), (96) 
C'il„+i=™"icS,„+iy. (97) 

^ (2) 

The z-component, ^ is found to be unnecessary here. From these results, we finally have 

-IC" (98) 



W = n\ 



f(0)A(2) , P(0)A(2) 



(b) Coupling through the lower path (n = n — 1) 

The expressions for C*r,+,n-i and C'^^+.„_i are the same as those given by equations (81) 
and (82), respectively, and Cl^]_ n-i is unnecessary, as in the case of (b) of subsubsection 7.1.1. 
Hence, all of the procedures go on in parallel, like in subsubsection 7.1.1. In subsubsection 
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7.1.1. in deriving from equation (49), [uj — {m + — nQ^ was replaced hj n"^ — nQ^, 

since the horizontal resonances are considered there. In the present case of vertical resonances, 
it can be approximated by n{kHyQ\/ (h^lj_ — k^). Hence, in the present case we have 



2n 



'^'+'"-^ a;-(m + l)l]^'+'"-^. 



X , ^ 7] + 



dr ' q:' 



drj 2Q rj 

dr uj — {m+ 1)Q r. 



(99) 



Ca,.-1 = ^(fci^)-^^%^CS,„-l. (101) 



^ ('2'! ^ f'2") 

Expressions for C^^l^n-i and 4. „_i are the same as those given by equations (91) and (92). 
Hence, we finally have 

W = n\ S,n-l(^S,n-l+(n-l)! gVl^iln-l = -35^ICS,n-ir. (102) 

7.2. P-Mode Oscillations 

In the case of p-mode oscillations, —{uj — niQ) + nQ']_ ~ 0, but the dispersion relation 
shows that it can be approximated as —n{kHyQ']_/[nQ']_ — k^). Hence, from equation (34) we 
have a relation between ^f'^ and ^IP\ which is 

^l''=^{kH)-^e^. (103) 

Furthermore, equation (33) gives 

i^;^=^^^6'^■ (104) 

These relations show that in the p-mode oscillations ' is major and that and are 
smaller than ^f'^ by a factor of kH, as expected. 

As mentioned before, we have only horizontal resonances, i.e., the vertical resonances 
are absent in the case of p-mode oscillations. 

7.2.1. Horizontal resonances of p-mode oscillations 

As in the case of g-mode oscillations, two coupling paths are treated separately, 
(a) Coupling through the upper path {h = n + 1) 

Expressions for C'^,+,n+i C'^]]_^n+i a^^^ the same as equations (79) and (80), respec- 
tively. {Cj^^_^_^n+i is unnecessary, since it is small.) The expressions for (^^l^^+i a^^d Ci^+,n+i 
also the same as those given by equations (83) and (85), respectively. The coupling terms. 
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C^J^^n_^_i are thus the same as those given by equations (86) and (87). The term 



.(2) 



I is unnecessary, since it is small and does not contribute to the final results. 



Using the above results, we finally have 



W 



(2) 



■,n+l 



(2) 



nl 



■IC 



(1) I 

r,+,n+l I 



(105) 



This expression for W is the same as equation (88), but its concrete figure is different from that 
of equation (88), since the expressions for and ^1°^ are different from those of the g-mode 
oscillations (see table 1). 

(b) Coupling through the lower path [n = n — l) 

The coupling terms in the first stage are found to be 



^(1) 



a 



,n-l 



2/(0) 



dr 



_ dn 
-rj + n— 
dr 



r 



(106) 
(107) 



These expressions for C^^^^^„_x and Cl^^j^^^-i different from equations (81) and (82), respec- 
tively, since the terms of ^f"* can be neglected in the present case. 
Using the above expressions, we have 



^(1) 



Expressions for „_i and Cz^+.n-i are not shown here, since they are the same as equations 
(84) and (85), respectively, except for the subscript n + 1 being replaced by n — 1. Consequently, 
as the final results concerning W , we have 



2^(0) J _ 



dr 



-rj + n 



dr] 
dr 



2Q 



V 



io — {m + l)f2 r 



(108) 



Ty = (n-l)!efC'f 



(2) 



(1) I 

r,+,n— 1 1 



(109) 



So far, we have restricted our attention to the coupling cases of m ^ m + 1 ^ m. Even 
in the case of m ^ m — 1 ^ m., however, the derivation of the W^s is similar. The results 



show that the final expressions for W^s are the same, except that the subscript + in 



or CXn-i ^ equatious (88), (93), (98), (102), (105), or (109) is replaced by 



:(i) 



(1) 

r,+,n+l 5 



:(i) 



the subscript — . The difference between (r^l^^+i ^-^d Cr,-,n+i5 fo^ example, is that m + 1 in the 
expressions for Cr,+,n+i is changed to m — 1 in (r,-,n+i- 

Next, we summarize the expressions of growth rate. As mentioned in section 5, the 
sign of [uj — {m±l)Q]/{uj — mQ) is negative in the horizontal resonances of g-mode oscillations. 
Hence, in this case we have, from equation (73) and expressions for W's given above. 



(1) 



-UJi 



-^2 



nfi2)2 



G 



H.±,n 



ijj — rriQ 



nl IC 



(1) I 



E 



(110) 



which is positive, and thus the resonant oscillations are amplified {—ui > 0). In the cases of 
horizontal resonances of p-mode oscillations (excluding n = 0), on the other hand, the sign of 
[u — {m ± 1)Q]/ {uj — mQ) at the resonant point is positive, and we have 
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E 



which is negative and the resonances lead to damping of the oscillations. 

In the case of vertical resonances of g-mode oscillations, the situations are a little dif- 
ferent. The sign of [cu — (m ± l)r2]/(c<j — mVt) is positive in this case. In the case of vertical 
resonances, the sign of J{l/D)dr is opposite to that in the case of horizontal resonances [see 



and the oscillations are damped. 
8. Discussion 

In this paper the resonant excitations of disk oscillations on relativistic disks deformed 
by a warp have been examined, considering non-linear couplings between the oscillations and 
the warp. The non-linear couplings that we consider are schematically shown in figure 1. The 
resonant processes involved in the couplings bring about amplification, or damping, of the 
original disk oscillations. First, a general stability criterion of the resonances is derived by 
using the Lagrangian formulation developed by Lynden-Bell and Ostriker (1974). Then, the 
criterion is applied to studying the resonant amplification. 

There are three types of resonances: i) horizontal and ii) vertical resonances of g-mode 
oscillations and iii) vertical resonances of p-mode oscillations. For each case, four coupling 
paths between the oscillations and the warp are possible, i.e., m — m ± 1 m and n n±l 
— s> n. Detailed examinations where the resonances occur and how much the frequencies of the 
resonant oscillations are were made in Paper I. In this paper we thus concentrate on whether 
these resonances amplify or dampen the oscillations. The results concerning their growth 
rates are summarized in equations (110)-(112). It should be emphasized that the presence 
of resonances is completely an effect of general relativity. This is because a difference in the 
spatial distributions of Q{r) and K,{r) is essential for the resonances, and it results from the 
disk being general relativistic. 

Although the detailed interaction processes are different in each type of resonances and 
in each coupling path, the final expressions for W have similar compact forms for all cases 
[see equations (88), (93), (98), (102), (105), and (109)], which suggests the correctness of our 
manipulations. [The results of Kato (2003) are not always correct, since there are some errors 
in the calculations.] 

In the case of horizontal resonances, W can be compactly expressed in terms of 
while it is expressed in terms of Q^^ in the case of vertical resonances, as expected. Although 



equations (69) and (72)]. The signs of the W's given above are also opposite 
case of horizontal resonances. Considering them, we have 



to those in the 



n/8 Gv,±,n {h-2)\\ci% 
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Table 1. Growth rates and resonant radii. 



Type of resonance 


Path 


Growth rate 


Radii 




g-mode, horizontal 


Ti. = n. -\~ 1 




g 




g-mode, vertical 


n = n+1 


—afl (damping) 


^ 3rg, 3.62rg 


6.46rg 




n = n — 1 




^ 3rg, 3.62rj 


6.46rg 


p-mode, horizontal 


n = n + 1 


—a{kH)'^n (damping) 


^ 3rg, 3.62rg 


6.46rg 


(n^O) 


n = n — 1, 


—aCl (damping) 


~ 3rg, 3.62rj 


r, 6.46rg 



the expressions for W are similar in all cases, concrete figures of the growth rates are different in 
each case, since the expressions for (or Q^^) are different in each case. Let us now estimate 
the order of the growth (or damping) rate. The amplitudes of intermediate oscillations, \ 
are given in section 7 for various types of resonances. Hence, using them we can estimate the 
order of the growth rates, which are summarized in table 1. In table 1, a denotes \T]\'^/rcL, 
where L is the radial size of the wave packet of the oscillations, and E has been estimated by 
E ~ l^^^^pL/r for g-mode oscillations and by ~ \Q^^\'^L/r for p-mode oscillations. The radii 
where the resonances for the oscillations of A; ~ occur are also given in table 1 (see also Paper 
I). It is noted that in table 1 p-mode oscillations of n = are included in g-mode oscillations. 

The radii given in table 1 are those where both curves specifying the propagation region 
and the resonance region cross on the u-r plane. In other words, they are radii where long- 
wavelength oscillations have resonances. We think that long-wavelength oscillations are of 
interest, since they are most observable (little phase mixing and little group velocity in the 
radial direction). The growth (or damping) rate of the oscillations, however, tends to zero in 
the limit of long wavelength of kH = 0. Hence, the oscillations that are of interest are those 
slightly apart from the radii given in the table. 

Table 1 shows that the vertical resonances of the g-mode oscillations with n = n + l and 
the horizontal resonances of the p-mode oscillations with n = n — l are strongly damped [see the 
absence of a small factor of (kH)]. The growing modes are horizontal resonances of the g-mode 
oscillations alone, including the horizontal resonances of the p-mode oscillations of n = 0. The 
growth rates of these oscillations are on the order of aQc times (kH)"^. This suggests that the 
quasi-periodic oscillations caused by our present mechanism are observable when the disks are 
moderately thick. 

Let us hereafter focus only on the growing oscillations, which are all amplified around 
4. Org. The reason why 4. Org is the resonance radius is that the resonance condition for horizontal 
resonances of g-mode oscillations is 2k = Q (see section 5), and this is realized there in the case 
of the Schwarzschild metric. The frequencies of the resonant oscillations are k, Q±k, 2Q±k,... 
at r = 4rg (see Paper I) , depending on the modes (m and n) of the oscillations. In units of 

° a in table 1 denotes \ri\l/rcL. 
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(GM/r^y^"^, these frequencies are 

0.0625, 0.1875, 0.313, 0.4375, (113) 

and their ratios are 1 : 3 : 5 : 7... The modes of oscillations realizing these frequencies are 
summarized in table 2. 

If we are concerned with what frequencies are observable in luminosity variations, some 
other frequencies should be added in the list of equation (113) for the following reason. Let us 
consider one-armed (m = 1) oscillations with frequency u. Then, a one-armed pattern revolves 
around the disk center with frequency u. After half a revolution, however, the arm appears on 
the another side of the disk. Hence, if we observe luminosity variations (not velocity variations) 
resulting from the whole disk, the phases of the luminosity variations become the same as the 
initial one after half a revolution. That is, we observe a luminosity variation with frequency 2u. 
As shown in table 2, the oscillations of k (or Q — n) are realized by one-armed oscillations. Hence, 
we can expect luminosity variations with frequency 2k(= Q). A similar argument shows that 
we can also expect a luminosity variation of 2{Q + k) or 2(2^2 — k), which is 6k,, since oscillations 
with Q + K, or 2Q — k, occur for one-armed oscillations. The above argument, however, cannot 
be applied to m = 2 oscillations. One-armed oscillations have a unique characteristic in this 
sense. Summing them, we can say that the observed frequencies are some hamonics of k, as 
follows: 

0.0625, 0.125, 0.1875, 0.313,.. (114) 

in the case of the Schwarzschild metric. Their ratios are 1: 2: 3: 5: 6: 7 ... 

The lowest frequency of oscillations, k, is hidden by the oscillations of 2k. The most 
observable pairs of oscillations are thus 2 : 3. This will be the reason why 2 : 3 oscillations are 
prominently observed as QPOs in black hole candidates. The above arguments suggest that 
other sets of oscillations, for example, 1 : 2 and 3 : 5, are also observable candidates. It is 
noted that the importance of 2 : 3 pair oscillations was first emphasized by Abramowicz and 
Kluzniak (2001) and Kluzniak and Abramowicz (2001), although they did not present definite 
resonant mechanisms. 

Next, we estimate the masses (and spins) of four black-hole objects with pairs of QPOs. 
In all of them, the pairs are close to 2 : 3. If we take the observed upper frequencies of these 
QPOs to be 3K(rmax) with the Schwarzschild metric, the masses of the central black hole can 
be calculated for each object. The results are shown on the fourth line of table 3. The second 
line of table 3 shows the observed pairs of frequencies in units of hertz, and the bottom line 
of the table shows the mass ranges (in units of the solar mass) which were evaluated from the 
observed data. In XTE 1550—564 and GRO 1655—40, the masses estimated by our model 
are slightly below the mass ranges obtained from observations. One of the reasons would be 
that effects of the rotation of the central objects are not taken into account in evaluating the 
frequencies. 
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Table 2. Frequencies and modes of resonant oscillations. 



Frequencies 




Modes of oscillations 






in units of {GM/ 


^g)^^^ g-mode (arbitrary 7i) p-mode (n = 0) 


K, Q 


— K 


0.0625 


m = 0, 1 


m = 0, 1 


n + i 




K U.iO/0 


m== 1, 2 


m = 1, 2 


2fi-|- 




K U.oio 


m = 2, 3 


m = 2, 3 






Table 3. Estimated i 


nasses of four black hole candidates. 


Objects 






XTE 1550-564(*) 


GRO 1655-40(*) H 1743-322(+' 


Frequencies 




^llO, luo j 


(184, 276) 


(300, 450) (160, 240) 


a 




IVl / IvIq 


M/Mq 


M/Mq M/Mq 


0.0 


4.0 


iz.u 


8.19 


4.75 8.90 


0.1 


3.84 


io.4 


8.78 


5.04 9.44 


0.2 


3.66 


14.3 


9.26 


5.37 10.1 


0.3 


3.48 


15.3 


9.93 


5.76 10.8 


0.4 


3.29 


16.5 


10.7 


6.22 11.7 


0.5 


3.08 


18.0 


11.7 


6.79 12.7 






10.0 - 18.0 


8.4 - 10.8 


6.0 - 6.6 not measured 



The processes for deriving the resonance condition 2k = f2 in this paper suggest that 
even when the central object has rotation, the condition for horizontal resonances of g-mode 
oscillations remains 2/t = f2 in the lowest order of approximations. Using the Kerr metric with 
the spin parameter a(0 < a < 1), we can calculate, as a function of a, the radius where 2n = Q 
is realized. Then, k (or Q) values at the radius are obtained as functions of the mass and spin 
of the central objects. Comparing the frequencies calculated in this way with observations, we 
can estimate the masses of the central objects as a function of a. The results are also presented 
in table 3, which shows that the masses evaluated by our model are in good agreement with 
observations for moderate values of a. Let us denote the upper frequency of the twin pairs by 
^upp- Then, the relation between i/upp and the mass, M, for a given a is 

^upp = A(^— j kHz, (115) 

where A = 2.1, 2.4, 2.8, and 3.4 for a = 0, 0.2, 0.4, and 0.6, respectively. The relation in the case 
of a = 0.4 is just the same as that obtained by McClintock and Remillard (2003) as a fitting 
formula of observations. 

The frequency ratios of some observed QPO pairs are slightly deviated from 2:3. This 
occurs natually in our model, since oscillations propagate in the radial direction as a wave 

° (*) McClintock & Remillard (2003), (+) Homan et al. (2003). See also Abramowicz et al. (2004b). 
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packet, changing the resonant frequencies and their wavelengths. Another possibihty of the 
deviation is precession of the warp. A precession of the warp is generally expected when the 
central object has rotation. If this interpretation is correct, the frequency of the precession 
is related to the deviation, and the precession might be observed as a low- frequency QPO. 
Large Q- values of QPOs are also expected in our model, since the oscillations are not trapped 
standing ones. Furthermore, different modes of oscillations with different paths of couplings 
appear with close frequencies, as shown in table 2 (see also Paper I). 

If our present model represents the real mechanism of kHz QPOs, observations of QPOs 
are a powerful tool for evaluating the mass and spin of the central black holes. Furthermore, 
it is implied that a warp occasionally appears in the innermost region of relativistic disks. 
This would be of importance for understanding the activities, structures, and evolutions of the 
innermost region of black-hole disks. 

Finally, it is noted that the basic process of our resonance model is applicable to any 
deformed disks. For example, if the disks are deformed by the spin of the central objects through 
a spin-disk coupling, we can expect disk oscillations whose frequencies are correlated to the spin 
frequency. It is known that in some X-ray binaries with neutron stars, the frequency difference 
between the high frequency QPO pairs is equal to, or half of, the spin frequency (Wijnands et 
al. 2003). In a subsequent paper we will examine whether the phenomenon could be explained 
by a simple extension of our present resonance model. 



Appendix 1. Expressions for C in Cylindrical Coordinates 

The coupling term C(^,^) consists of two terms as C = C^, + Cp. Expressions for 
and Cp in cylindrical coordinates are 



where 



PoCiP^ip 



-^Po 



Q^2 



drdz 



dz^ / dr dr r dr 



2 



-po 



/ dr\r dr J r drdz 



PoC^,z = --po 



PoCp^r - 

PoCp^ip 

PoCp^z 



dT„ 



,d^ 

dT^r 



d^ 

drdz 
dTzr 1 



d"^ \ dipo 



dr rdip 
dTr^ dT, 



dz 



T 



-T. 



T 



Po 



dr 
dTrz 
dr 

dr 



rdip 



dz r 



■T. 



ipr Ji 



dT^z dTz, 



rdip 



\rd(p 



dz 



+ 



d^r diz 



dr dz dr 



^d'lpo 



dz'^ J dz r drdz 



(116) 
(117) 
(118) 
(119) 
(120) 
(121) 

(122) 
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T 



T 

J. r 



T. 



ipr 



T. 



X 



T 



- Zip 



and 



^ 7. 



Po 



Po 



Po 



Po 



Po 



Po 



Po 



dr rdip r 



rdif 



r 



+ 



dz rdif 



dr dz \rdip 

dr dr \ rdip 
'dip f djr 



di, di. 



r J dz 



+ 



dz dz 
dip, di. 



r J dr dz dr 
r 



rdif 



dz rdif 



. dr dz 



rdif 



r J dz dz dz 



r ^d^d^^d^d^ 
dz dr 



dr dr rdip dr 
dr \rdip r J rdip \rdip r J 



diz di. 



dz rdip 



diz dir , di^ di^ , di^ di 



+ 



dr dz rdif dz dz dz 



Appendix 2. Expressions for C'^.^ n' ^^j^i.n' ^'tX,n 
From the expressions in appendix 1, we have 



^r,+,n+l-^l^^[ 

d^r dr] , d^^ rj dlnH 
dr dr dr r dr 



1 rf2(]2 



+n{n — 1) 



dr r 



n———-irV ] T7^ 

d\nH 



n d 



dr H rpooH dr 



rpooir 



dr 



dr r 



^. i^ 
+ in— 

r 



fm + 1)- 

dr r 



dr 
drj 



d 



rdr ^ ^ J.2 dr ^ r ^ 



■n — : — ip>- + i 



n 



dr 



r PooH^ ^ dr 



dlnH - r] 

r 

r 

PooH'^-Xv 



^ rj 

dlnfi^ - _ H - / ^ drj rj 



drj 
' dr 



'i{n — l) 



Hr 



^(1) _ r)2 

^z,+,n+l — i'i 



dr 
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The above expressions for the C"s are those for the path ofm^m+l-^m (notice + in the 
subscript). The expressions for the C"s resuhing from the path ofm— s>m— 1 ^ m (i.e, C"s 
with the subscript of — ) are similar to those with a + sign. That is, expressions for C's with 
the - sign are obtained by changing i and m in the expressions for C's with the + sign to —i 
and — m, respectively. Furthermore, rj is changed to 77*. 

Appendix 3. Expressions for C'^,±,n? ^!pX,ni ^i%,h 

^ (2) 

Expressions for C^±n can be derived from G's given in appendix 1 by similar procedures 
as in appendix 2. The results of manipulations show that C*^ + ,^+i(C+,li+i! ''7*) is related to 
C'r,-,n-i(C-,n-i)''7*)- That is, the former expression is obtained from the latter one by changing 
the argument, in the latter to C+L+i ^i^d futher m and n in the latter, respectively, to 

m + 1 and n + 1. Similarly, C'r|+,n-i(C+,L-i!''7*) is obtained from C^r^_^n+i{$,^-]n+i-,'n*) by changing 
$!~-,n+i ill the latter to ^^'^ ^"^^ n to m+ 1 and n — 1 , respectively. 

The expression for C^}_^n+i{C^-]n+i^'n) is obtained from C*r,+,n-i(^+,n-i) ■'?) by changing 
^+.n-i to C-.li+i ^-iid m and to m — 1 and ra + l, respectively. Similarly, C*r,-,n-i(C-,L-i)''7) is 
obtained from C^]^^n+i{^^+,n+ii'n) by changing to C-,n-i m and n, respectively, to 

m — 1 and n — 1 . 

Expressions for C*^^j_^^ and C'^^].^^ are neglected here, since they are easily supposed from 
the above procedures. 



27 



References 



Abramowicz, M. A. & Kluzniak, W. 2001, A&A, 374, L19 
Abramowicz, M. A., & Kluzniak, W. 2003, astro-ph/0312396 

Abramowicz, M. A., Kluzniak, W., McClintock, J. E., & Remillard, R. A. 2004a, astro-ph/0402012 
Abramowicz, M. A., Kluzniak, W., Stuchlik, Z., & Torok, G. 2004b, astro-pli/0401464 
Hirose, M., k Osaki, Y. 1990, PASJ, 42, 135 

Homan, J., Miller, J. M., Wijnands, R., Steeghs, D., Belloni, T., van der Klis, M., & Lewin, W. H. G. 

2003, Atel 162 (http://integral.rssi.ru/atelmirror) 
Kato, S. 2001, PASJ, 53, 1 
Kato, S. 2003, PASJ, 55, 801 
Kato, S. 2004, PASJ, 56, 559 (Paper I) 

Kato, S., Fukue, J., & Mineshige, S. 1998, Black-Hole Accretion Disks (Kyoto: Kyoto Univ. Press) 
Kato, S., & Tosa, M. 1994, PASJ, 46, 559 
Kato, S., & Unno, W. 1967, PASJ, 19, 1 

Kluzniak, W., & Abramowicz, M.A. 2001, Acta Phys. Pol. B32, 3605 
Lubow, S. H. 1991, ApJ, 381, 259 

Lynden-Bell, D., & Ostriker, J. P. 1967, MNRAS, 136, 293 
McClintock, J. E., & Remillard, R. A. 2003, astro-pli/0306213 
Okazaki, A. T., Kato, S., & Fukue, J. 1987, PASJ, 39, 457 
Paczynski, B., & Wiita, P. J. 1980, A&A, 88, 23 
Tosa, M. 1994, ApJ, 426, L81 

van der Klis, M. 2000, Ann. Rev. Astron. Astrophys. 38, 717 
Whitehurst, R. 1988, MNRAS, 232, 35 

Wijnands, R., van der Klis, M., Homan, J., Chakrabarty, D., Markwardt, C. B., &: Morgan, E. H. 
2003, Nature, 424, 44 



28 



